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We  consider  the  stirring  of  an  inviscid  fluid  caused  by  the  locomotion  of  bodies  through  it.  The 
swimmers  are  approximated  by  non-interacting  cylinders  or  spheres  moving  steadily  along  straight 
lines.  We  find  the  displacement  of  fluid  particles  caused  by  the  nearby  passage  of  a  swimmer  as 
a  function  of  an  impact  parameter.  We  use  this  to  compute  the  effective  diffusion  coefficient  from 
the  random  walk  of  a  fluid  particle  under  the  influence  of  a  distribution  of  swimming  bodies.  We 
compare  with  the  results  of  simulations.  For  typical  sizes,  densities  and  swimming  velocities  of  krill, 
the  effective  diffusivity  in  this  model  is  five  times  the  thermal  diffusivity.  However,  we  estimate  that 
viscosity  increases  this  value  by  two  orders  of  magnitude. 


Munk  [1]  was  the  first  to  ask  whether  biology  has  an 
important  impact  on  mixing  in  the  oceans.  Since  mix¬ 
ing  affects  the  global  circulation  and  stratification  of  the 
oceans,  it  is  of  great  interest  to  physical  oceanographers 
to  settle  this  question.  Dewar  et  al.  [2]  have  proposed 
that  the  mechanical  energy  delivered  by  the  swimming 
motions  of  the  marine  biosphere  could  amount  to  almost 
1012  W,  a  figure  comparable  to  the  energy  delivered  by 
the  winds  and  tides.  This  suggests  a  biological  origin  of 
about  33%  of  the  mixing  in  the  oceans,  an  enormous  fig¬ 
ure.  By  assuming  that  this  energy  is  delivered  to  the  top 
three  kilometers  of  the  oceans,  they  estimate  an  effective 
diffusivity  produced  by  swimmers  to  be  approximately 
.2  cm2/  sec,  or  about  100  times  the  molecular  value  for 
heat.  Kunze  et  al.  [3]  have  measured  elevated  levels  of 
ocean  turbulence  due  to  swimming  krill,  though  some 
have  questioned  whether  this  turbulence  can  efficiently 
overturn  a  stratified  medium  [4,  5].  Katija  and  Dabiri  [6] 
suggest  that  the  displacement  of  fluid  particles  by  swim¬ 
ming  bodies,  which  viscous  effects  can  lengthen,  is  more 
relevant  to  stirring  than  the  scale  of  turbulence  they  pro¬ 
duce.  This  is  the  viewpoint  adopted  in  this  Letter. 

Huntley  and  Zhou  [7]  considered  the  energy  produced 
by  11  representative  species  of  schooling  animals,  from 
krill  to  whales,  and  found  that  regardless  of  size  the  en¬ 
ergy  input  to  the  ocean  per  unit  mass  from  the  swim¬ 
ming  of  these  animals  was  roughly  a  constant  of  order 
10-5  W /  kg.  If  the  average  the  biomass  of  the  oceans  has 
a  volume  density  of  c  over  a  volume  V  in  cubic  meters,  we 
can  thus  (assuming  biological  materials  have  the  density 
of  water)  arrive  at  a  total  energy  input  of  10~2cV  W. 
The  area  of  the  oceans  is  about  3.6  x  10 14  m,  and  if  we 
assume  the  biomass  is  evenly  distributed  to  a  depth  D 
meters,  we  get  an  input  of  3.6Dcx  1012  W.  This  suggests 
that  for  a  typical  depth  of  one  kilometer  the  volume  den¬ 
sity  of  the  biomass  should  be  of  order  10_3-10-4,  and 
thus  that  the  organisms  form  a  dilute  suspension,  even  if 
their  distribution  is  patchy  and  organized  into  schools. 

The  object  of  the  present  paper  is  to  determine  the  ef¬ 
fective  diffusivity  and  the  statistics  of  the  concentration 


field  of  a  passive  scalar  that  can  result  from  the  fluid  mo¬ 
tions  caused  by  such  dilute  arrangements  of  swimming 
animals.  The  scalar  could  be  for  example  heat,  salt,  or 
nutrients,  and  is  assumed  to  have  negligible  feedback  on 
the  flow  (at  least  at  small  scales) .  The  animals  considered 
by  Dewar  et  al.  [2]  and  Huntley  and  Zhou  [7]  have  large 
Reynolds  number,  typically  102-107.  Our  focus  is  there¬ 
fore  distinct  from  the  mixing  that  can  occur  from  dense, 
tightly  interacting  suspensions  of  Stokesian  swimmers, 
as  in  [8],  though  it  shares  many  common  features  such 
as  linearity  in  the  density  of  swimmers.  Real  swimmers 
presents  a  wide  variety  of  motions  and  wakes,  but  these 
have  generic  forms  in  the  far  field  [9].  For  example,  the 
far  potential  field  of  a  neutrally  buoyant  fish  in  steady  un¬ 
accelerated  swimming  decays  at  least  like  a  quadrapole. 
Following  [6],  we  will  model  the  swimmers  by  identical 
cylinders  or  spheres  moving  in  potential  flow,  but  these 
are  merely  stand-in  examples.  We  emphasize  that  though 
the  motivating  application  comes  from  oceanography,  the 
simple  model  we  introduce  can  be  applied  to  a  range  of 
systems,  such  as  the  mixing  caused  by  vortices  or  sus¬ 
pensions  of  solid  particles. 

Dilute  suspensions  of  swimmers.  As  an  extremely  sim¬ 
ple  model  of  stirring  by  swimmers,  we  consider  a  fluid 
particle,  called  the  target  particle,  which  is  influenced  by 
the  occasional  passing  of  swimming  bodies.  We  assume 
that  every  swimmer  moves  at  a  fixed  speed  U,  over  dis¬ 
tances  O (L)  large  compared  to  the  “range  of  influence.” 
This  range  will  typically  be  a  few  body  lengths  normal 
to  the  swimming  path,  and  represents  the  distance  from 
a  target  particle  where  interaction  with  the  swimmer  be¬ 
comes  significant.  We  assume  that  encounters  of  our  tar¬ 
get  particle  with  swimmers  are  occasional,  in  that  at  each 
encounter  the  particle  moves  by  a  distance  A  (a),  where 
the  impact  parameter  a  is  the  perpendicular  distance  of 
the  initially  unperturbed  particle  from  the  extended  line 
of  motion  of  the  approaching  swimmer.  Since  we  are  as¬ 
suming  that  L  is  large  compared  to  the  distance  of  signif¬ 
icant  interaction,  we  may  assume  that  by  an  “encounter” 
we  mean  that  the  A  (a)  may  be  computed  from  motion 
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of  the  swimmer  along  a  doubly-infinite  line.  Each  en¬ 
counter,  the  fcth  say,  moves  the  particle  a  distance  A  (a*,), 
in  the  direction  of  the  unit  vector  fk.  Thus  after  M  en¬ 
counters  the  position  of  a  target  particle  initially  at  x0 
is  x m  given  by 

M 

XM  =  X o  +  ^  a {ak)rk  .  (1) 

k= 1 

Assuming  infrequent  encounters,  we  regard  the  ak ,  rk 
as  independent  and  identically  distributed  random  vari¬ 
ables.  Then  [10] 

(|a;|2)  =  M( A2(a))  +  (J2  A(aj)A(ak)  rj  ■  rfc),  (2) 

where  we  assume  spatial  homogeneity  to  set  Xq  =  0,  and 
the  angle  brackets  denote  ensemble  averaging.  The  fk 
are  isotropically  distributed,  so  the  second  term  above 
vanishes. 

To  evaluate  M(A2(a)),  consider  a  swimmer  at  a  large 
distance  L  from  the  target  particle  at  the  origin  O. 
Imagine  a  “target  disk”  (or  sphere  in  3D)  of  radius  R, 
with  I  <  J?.  C  1,  where  t  is  a  typical  length  scale  for  the 
swimmer.  The  fraction  of  swimmers  that  will  hit  the  tar¬ 
get  disk  (sphere)  is  2R/2nL  (ttR2 /4ttL2  in  3D).  We  can 
use  the  number  density  and  the  volume  of  the  ‘shell’  at  a 
distance  L  to  find  that  the  number  of  swimmers  hitting 
the  target  from  a  distance  L  is  2 RncLL  ( TrR2ndL  in  3D), 
where  n  is  the  number  density  of  swimmers.  We  now 
integrate  from  0  to  LT,  since  swimmers  further  than  Ut 
cannot  hit  the  target,  to  find  the  number  of  swimmers  M 
that  will  hit  the  interaction  disk  (sphere)  in  time  t:  we 
find  M  =  2RUnt  ( M  =  ■nR2Unt  in  3D).  The  expression 
for  the  squared  displacement  (2)  is  now 
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FIG.  1:  Displacement  A  (a)  due  to  a  cylinder,  as  a  function  of 
the  impact  parameter  a.  The  dashed  lines  are  the  asymptotic 
forms  for  small  and  large  a  of  Eq.  (5).  A  few  typical  trajecto¬ 
ries  are  overlaid:  for  large  a  the  trajectory  is  almost  circular, 
but  is  still  ribbon-shaped. 

for  t  ranging  from  — oo  to  oo.  The  impact  parameter, 
a,  appears  as  the  initial  vertical  position  of  the  target 
particle.  This  classical  ‘drift’  problem  has  been  treated 
in  great  detail  by  many  authors;  see  for  example  [11,  12]. 

As  a  working  example  of  the  displacement  due  to  a 
moving  body,  we  shall  first  treat  the  cylinder.  We  con¬ 
sider  the  displacement  in  two-dimensional  flow  due  to  the 
passage  of  a  cylinder  of  radius  t  moving  at  speed  U  in 
potential  flow.  Equation  (4)  can  be  integrated  numeri¬ 
cally  for  a  given  a  to  obtain  A(a),  plotted  in  Fig.  1.  The 
characteristic  ‘ribbon’  shape  is  evident.  Also  shown  are 
the  asymptotic  forms  [11] 

4(«)/<  +  “<<f;  (5) 

[7r(£/a)3/2,  a>f. 


rR 

(|a;|2)  =  2M  /  p(a)A2(a)da 
Jo 

since  the  largest  value  of  the  impact  parameter  is  R. 
Here  p(a)  is  the  probability  distribution  of  impact  pa¬ 
rameters;  since  the  swimmers  are  assumed  to  arrive  from 
far  away,  this  is  p(a)  =  1/R  ( p(a )  =  2na/nR2  in  3D). 
Combining  these  results  and  taking  R  — >  oo,  we  find 

2  (2UntJ0xA2(a)da  (2D); 

Vl  '  '  ylnUnt  a  A2 (a)da  (3D),  1  j 

assuming  that  the  integrals  converge.  Since  the  effective 
diffusivity  k  is  defined  by  ( | aj | 2)  =  2 dut,  with  d  the  spatial 
dimension,  Eq.  (3)  can  be  used  to  determine  n. 

Displacement  due  to  a  moving  body.  To  find  the  dis¬ 
placement  A  (a)  of  a  target  particle  due  to  a  swimming 
body  coming  from  infinitely  far  away  and  swimming  in 
the  z  direction,  we  need  to  integrate  the  equations  for 
the  position  of  a  target  particle  at  x(t), 

x  =  u(x,t) ,  x(— oo)  =  (a,  0,  —  oo)  (4) 


There  is  an  integrable  logarithmic  singularity  for  small  a. 
For  large  a,  the  trajectories  are  almost  circles.  By  com¬ 
bining  these  asymptotic  limits  with  numerical  integra¬ 
tion,  we  find  the  integral  from  Eq.  (3)  is  J0°°  A 2  (a)  da  ~ 

2.371'3.  But  note  that  A2(a)da  ~  2.31£3:  the  integral 
is  completely  dominated  by  “head-on”  collisions  (97%  of 
the  integral). 

For  a  sphere  in  three-dimensional  potential  flow,  the 
displacement  A  (a)  appears  superficially  much  as  for  the 
cylinder  in  Fig.  1,  but  drops  off  more  rapidly  for  large  a: 


A {a) /l 


|  log(£/a)  —  .582,  a  <C  l\ 

97r(£/a)5/64,  a><. 


(6) 


As  for  the  cylinder,  there  is  an  integrable  logarithmic 
singularity  for  small  a,  and  the  integral  in  Eq.  (3)  is 
fQ  aA2(a)da  ~  .254£4.  The  small  a  logarithmic  singu¬ 
larity  is  mollified  by  the  extra  factor  of  a  in  the  integral, 
but  still  aA2(a)da  ~  .250£4:  as  it  was  for  for  the  cylin¬ 
der,  the  integral  is  completely  dominated  by  “head-on” 
collisions  (98%  of  the  integral),  due  to  the  rapid  decay  of 
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FIG.  2:  (a)  Initial  configuration  of  N  =  100  cylinders  moving 
at  constant  speed  U  =  1  in  random  directions.  Each  cylinder 
has  radius  1=1  and  the  periodic  box  size  used  in  the  simu¬ 
lation  is  L  =  1000  (cylinders  not  to  scale),  (b)  Trajectory  of 
target  particle,  initially  at  the  origin,  integrated  for  10°  time 
units. 

the  displacement  with  impact  parameter.  For  both  the 
cylinders  and  spheres,  the  logarithmic  singularity  is  the 
dominant  contribution  to  the  integral.  The  coefficient 
of  the  logarithm  in  Eqs.  (5)  and  (6)  is  given  by  the  lin¬ 
earized  flow  near  the  stagnation  points  at  the  front  and 
rear  of  the  cylinder  or  sphere,  suggesting  that  the  integral 
is  easy  to  approximate  for  more  complicated  swimmers. 
Putting  all  the  numerical  factors  together,  we  find  the 
effective  diffusivity 

f  1.19  Un£3  (cylinders); 

.266  Uni4  (spheres). 

Direct  simulation  of  dilute  suspensions.  To  validate 
our  theoretical  predictions,  we  consider  the  encounters  of 
2D  swimmers  moving  in  straight  lines  within  the  square 
S:  —\L  <  x,y  <  \L.  We  assume  each  swimmer  is 
a  cylinder  of  radius  £  =  1,  with  L  £.  We  initially 
place  N  swimmers  at  random  positions  within  S ,  which 
subsequently  move  with  unit  speed  in  a  random  direc¬ 
tion.  (Figure  2(a)  shows  a  typical  initial  configuration.) 
Positions  are  subsequently  computed  mod  L  in  both  di¬ 
rections,  maintaining  the  number  density  n  in  S.  Dilute¬ 
ness  requires  that  nt 2  -C  1.  A  target  particle  initially  at 
the  origin  moves  under  the  potential  flow  created  by  all 
of  the  cylinders.  Since  the  cylinders  are  typically  well- 
separated,  we  compute  their  net  velocity  field  by  linear 
superposition.  We  show  an  example  of  this  computa¬ 
tion  in  figure  2(b),  for  the  values  L  =  1000,  n  =  10-4. 
The  larger  “ribbons”  caused  by  drift  are  easily  identified, 
suggesting  that  in  this  dilute  limit  the  approximation  of 
encounters  as  being  independent  will  hold. 

Figure  3  shows  the  mean-squared  displacement  (|a:|2) 
of  a  target  particle  over  2  x  106  trials  (realizations) 
for  N  =  10  cylinders,  again  with  unit  radius  and  speed. 
The  solid  line  confirms  that  ( |  a;  | 2)  grows  linearly  with 
time,  and  the  dashed  line  shows  the  2D  theoretical  pre¬ 
diction  (7)  for  cylinders.  The  discrepancy  is  due  to 
Eq.  (7)  only  being  valid  in  the  limit  of  infinite  dilution, 
i.e.  as  n£2  — >  0.  We  have  verified  that  as  the  dilution 
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FIG.  3:  The  mean-squared  displacement  (solid  line)  of  a  tar¬ 
get  particle  for  2  x  10G  realizations  of  IV  =  10  cylinders,  with 
otherwise  the  same  parameters  as  in  the  caption  to  Fig.  2(a). 
The  dashed  line  shows  the  squared-displacement  predicted  by 
Eq.  (3)  and  (7),  using  a  number  density  n  =  10-5.  The  dis¬ 
crepancy  between  the  two  lines  is  due  to  (7)  being  valid  in 
the  limit  of  infinite  dilution.  The  background  shows  a  few 
individual  trajectories:  the  mean  is  dominated  by  rare  large 
events. 

is  increased  the  slope  approaches  the  theoretical  predic¬ 
tion.  In  the  background  of  figure  3  we  show  the  squared 
displacement  for  some  individual  trajectories:  it  is  clear 
that  the  mean  is  dominated  by  a  few  trajectories  with 
large  displacements,  corresponding  to  small  impact  pa¬ 
rameter  a.  Analogous  remarks  apply  to  a  suspension  of 
spheres  in  three  dimensions:  the  theoretical  predictions 
are  verified  there  as  well. 

Typical  numerical  values.  We  will  use  values  for  typ¬ 
ical  krill  as  in  [4].  We  consider  spheres  of  radius  £  = 
1  cm,  swimming  speed  U  =  5  cm/ sec,  and  number  den¬ 
sity  n  =  5  x  10-3  cm-3.  Equation  (7)  then  gives  an 
effective  diffusivity  of  7  x  ICG3  cm2/ sec,  about  five  times 
the  thermal  molecular  value  1.5  x  10-3  cm2/  sec,  and  five 
hundred  times  the  molecular  value  1.6  x  10-5  cm2/ sec 
for  salt.  This  implies  a  considerable  enhancement  to  the 
molecular  diffusion,  but  we  emphasize  that  these  values 
apply  within  a  school  of  krill:  the  distribution  and  size  of 
the  schools  themselves  is  a  more  complicated  matter  [7]. 
Note  also  that  a  small  change  in  the  swimmer  size  has  a 
huge  impact:  for  a  radius  of  .5  cm,  the  effective  diffusiv¬ 
ity  is  4  x  10-4  cm2/ sec,  an  order  of  magnitude  smaller 
than  for  1  cm. 

Effect  of  viscosity.  We  expect  viscosity  to  greatly  en¬ 
hance  k.  This  will  be  the  focus  of  future  investigation, 
but  for  now  we  present  a  rough  estimate  of  the  impact 
of  viscous  no-slip  boundary  conditions  at  the  surface  of 
the  swimmer.  For  inviscid  flow,  the  displacement  func¬ 
tion  A  (a)  has  a  logarithmic  singularity  near  the  axis  of 
swimming;  for  viscous  flow  near  a  no-slip  boundary,  the 
displacement  function  has  the  stronger  singularity 

A  (a)  ~  Cl2 /a,  a  <  £,  (8) 

where  C  =  a/2/3  tt  for  a  sphere  of  radius  i  [13].  This  im¬ 
plies  that  the  3D  squared-displacement  integral  in  (3)  di- 
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verges  as  a  — >  0.  This  unphysical  situation  can  be  reme¬ 
died  by  capping  the  displacement  function  at  a  maximum 
value  A,  which  is  the  path  length  of  swimming.  In  other 
words,  a  particle  which  is  directly  in  the  path  of  the  swim¬ 
mer  cannot  travel  further  than  the  swimmer  itself:  this 
regularizes  the  integral  (3)  to  give 

pOO 

{\x\2)  ~  2nUnt  /  aA2(a)da  (9) 

J  Cl2  /  A 

where  the  lower  bound  of  the  integral  is  the  value  at 
which  the  displacement  (8)  achieves  its  maximum  allow¬ 
able  value,  A.  Assuming  that  most  of  the  contribution 
to  the  integral  arises  from  small  a,  as  was  the  case  for 
potential  flow,  we  have 

f  Uni4C2log(\/i).  (10) 

For  spheres,  this  is  k  ~  (2tt3 /9)Un£4  log(X/£)  = 
6.89  Uni4  log(A/£).  Inserting  the  same  numerical  val¬ 
ues  for  krill  as  before,  with  a  path  length  A  =  100  cm, 
we  find  k  ~  .8  cm2/ sec,  about  500  times  the  molecular 
value.  Thus,  including  the  effect  of  viscosity  and  finite 
path  length  has  increased  the  effective  diffusivity  by  a  fac¬ 
tor  of  100  over  the  inviscid  flow  case.  We  emphasize  that 
this  is  a  rough  estimate.  The  path  length  (or  swimming 
correlation  length)  is  a  measure  of  how  much  a  swimmer 
tends  to  move  in  one  direction  before  turning.  Our  cho¬ 
sen  value  of  1  m  is  not  based  on  any  evidence,  but  n  has 
only  a  weak  logarithmic  dependence  on  A.  Assuming  a 
value  A  =  10  m  only  raises  k  from  .8  to  1.2  cm2/ sec. 

Conditions  where  the  number  density  is  higher,  such 
as  in  aquaculture  [14]  or  locally  in  the  ocean  [3],  should 
also  increase  the  effective  diffusivity  well  above  molecu¬ 
lar  values.  The  simplicity  of  our  model  means  that  pref¬ 
actors  and  scalings  can  be  computed  accurately.  The 
numerical  constants  we  obtained  depend  mostly  on  the 
flow  near  the  stagnation  points  around  the  swimming 
body.  There  are  many  effects  that  we  neglected  but  must 
be  incorporated  in  a  more  complete  theory:  viscosity, 
boundary  layers,  more  realistic  shape  distributions  for 
the  swimming  bodies,  wakes,  spatial  correlations  between 
the  swimmers,  patchiness  and  schooling,  finite  correla¬ 
tion  length  of  swimming,  distribution  of  velocities  and, 
most  importantly  for  oceanic  applications,  buoyancy  and 
stratification  effects.  This  last  item  is  probably  the  most 
important:  stratification  can  cause  fluid  parcels  to  return 
to  their  initial  height  after  being  displaced  if  they  can’t 
equilibrate  their  density  with  their  surroundings.  Our 


simple  model  can  serve  as  a  platform  on  which  to  build 
complexity,  or  could  be  applied  to  other  fluid-dynamical 
systems  where  a  collection  of  objects  causes  mixing. 
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